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ABSTRACT 


This thesis describes the eonve zion of a computer 
program from Fortran IV for the NOS 1.2 operating system of 
the CYBER 175 cr CDC 6600 computer to Fortran IV compatible 
with the Naval Pcstgraduate School IBM 3033 system. The 
converted program, called FLO27, calculates the inviscid, 
three-dimensicnal transonic potential flow over wings or 
witng-body ccmbinations. The data input to FLO27 is 
extensive; therefore, an interactive program was developed 


to aid the user in building the required input data file. 
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I. INTRODUCTION 


In the Aeronautical Engineering curriculum graduate 
level aerodynamics course, AE-4501, the students are exposed 
to two computer crograms. One of these, prepared by the 
Douglas Aircraft Company, analyzes the potential flow around 
three-dimensicnal wings but is limited to incompressible 
eLow (Ref. 1]. The other progran, prepared by Cebeci, 
@atculates .the. friction drag for two dimensional 
incompressible flow over airfoils [Ref. 2]. A serious 
defect of these programs is that they are not 
state-of-the-art computer frograms. The Douglas program 
does not consider th2 effects of compressibility and the 
boundary layer program, in additicn to being restricted to 
incompressible flow, does not predict the laminar to 


turbulent transition point. 


A. EBACKGRCUND 

In 1980 the Department of Aeronautics at the Naval 
Postgraduate School acquired an intricat? computer program 
recently develoved by the Bo?ing Commercial Airplane 
Company. This state-of-the-art program calculates 


three-dimensional transonic flow over wings and hodies in 


both «he cuzer-inviscid flow regicn governed by «he 
transonic vrotential equation and the thin layer in which th2 
£2 rst.order, ccmpressible boundary layer equations are 
assumed tc te valid. 

The Boeing »vrogram as received was designed to be 
executed on a CDC 6600 or a CYBER 175 computer and was 
written using CDC FORTRAN IV extended language. This thesis 
therefore was primarily ccncerned with the conversion of the 
program to FORTRAN IV extended ccmpatible with the Naval 
Postgraduate School's (NPS) IBM 3033 system. The large 
modular orcgram was divided so that th2 potential flow 
analysis scrtien could be run separately. Simplified 


instructions for use of the program were also prepared. 


B. VISCOUS/INVISCID SYSTEF OF PROGRAMS 

The Viscous/Inviscid Wing System (VIWS) of programs 
calculates three-dimensional transonic flow over wings and 
wing kody ccmbinations including details of the laminar or 
turbulent flow in th2 three-dimensional viscous boundary 
layer. The flew field is calculated in two overlapping 
regions: an outer inviscid flow region governed by the 
transonic fotential equation, and a thin boundary layer in 


which the first order, three-dimensional, compressible 
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boundary laver equaticns are assumed to hold and in which 


+he effects of surface teat and mass transfer can be 
computed. A list of the VIWS of programs is presented in 


table I. 


TABLE I 


Viscous/Inviscid Wing Syst2m of Programs 


F1027 Jameson-Caughey inviscid, 
transonic wing code 

A411IN Reads gecmetry & velocity 
data, constructs coordinate 
systen 

VWIN Pcetential flew boundary 
layer interface 

A411AC1 Three-dimensional boundary 
layer pregran 

INTERP eee | layer potential 
flow interface 


Graphics display programs 


The basic sequence of calculations used by the VIWS to 
obtain matched viscous and inviscid solutions consists of an 
iterative looc in which the inviscid outer flow analysis and 
the beundary layer analysis are performed sequentially. The 
iterative sequence is ccntinued until either convergence 
(satisfactory matching) is achieved, or the maximum number 
of iterations sp¢écified by the user has been performed. The 
VIWS programing sequence is shown schematically in Fig. 1.1. 
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Final boundary 
layer calc? 


Generate boundary layer grid 
1st interface program 


Calculate final boundary layer flow 
1. Boundary layer program—upper surf. 
2. Boundary layer program—lower surf. 


Calculate transonic potential flow 
FLO 27 
Generate boundary layer grid 
Ist interface program 


Calculate 3-D boundary layer flow 
1. Boundary layer program—upper surf. 
2. Boundary layer program—lower surf. 


Plotting programs 


Set up 5” for next iteration 
2nd interface program 

1. Interpolation 

2. 5” under-relaxation 
3. Convergence test 


5° Yes 
converged 


Figure 1.1. Viscous/Inviscid Interaction Procedure 


The potential flow is calculated for the bare wing during 
the first iteration. In subsequent iterations, the effect 


of the boundary layer flew on the outer inviscid flow is 


wttvant 


worms GE 


‘ee ye 


felt as a modification to the wing shape through the 
addition of the boundary layer displacement thickness. 
Convergence is recognized, and the iterations are stopped, 
when the maximum change between the new and old displacement 
thickness, expressed as a fracticn of the maximum 
displacement thickaess, is less than the convergence 
tolerance chcsen by the user. 

The VIWS utilizes the Jameson-Caughey transonic inviscid 
mang program FLO2Z7,. to carry cut the potential flow 
analysis. The boundary layer analysis is performed by a 
finite difference boundary layer prediction program 
developed by the Boeing Commercial Airplane Company. The 
basic theory rehind the bcundary layer program is contained 
in [{Ref. 3]. A detailed descripticn of tha VIWS of programs 
(2xcluding the potential flow program FLO27) is contained 
in (Ref. 41. A basic guide to the use of the VIWS of 


programs is contained in [Ref. 5]. 
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II. POTENTIAL FLOW PROGRAM FLO27 


Eecause of the extensive length and number of program 
modules in the VIWS, the Fotential Flow Program, FLO27, was 
singled out for ccnversion. It was anticipated that FLO27 
would be run separately at first and recombined with the 
other crogram modules at scme later date when these modules 
were themselves translated for execution on the IBM 3033 


computer. 


A. RE-PROGRAMING 

The Potential Flow Pregram, hereafter called FLO27, was 
received cn magnetic tape and loaded into «he IBM 3033 mass 
stcrage system using the Jcb Control Languéage@ (JCL) routines 
presented in Arpendix A. The magnetic tape contained twenty 
m0) total files in’ which the format was 9 track, 1600 CPI, 
unlakled. The card image format for the sixteen (16) 
program files is 80 characters per record and the four (4) 
Output files contain 150 characters per record. The program 
amd output files on the criginal CDC tape are listed in 
fabie II. 

The FLO27 program was converted to FORTRAN IV extended 


Suitable for execution on the IBM computer using the NPS CDC 


14 


5 foe dttpnel- 


volt. Setervste® 


Pee 


Lelyeese® 


TABLE II 


CDC Magnetic Tape Files 


72550 FLO27 Potential Flow —-- 


Zz /3194 A4Y11IN Reads gecmetry velocity data 
constracts coor anaes systéen 

a7 3718 VWIN Potential flow boundary-layer 
interface 

4 /6611 A411A01 ineee -dimensional boundary 
layer program 

SO 719A INTERP CLAS Gs aioe potential flew 
interface 


Boeing McLean computer program 
[a6 700 | staepua| stzeanwise plots SSC—*d 


OUT4YI1L | Output from boundary-layer, 
lower surface 
OUT411U | Output from poundary-layer, 
upper surface 


eo rem conversicn guide [Rfef. 6]. The first step taken 


Procedure files 


consisted cf program compilation using the WATFIV compiler 


with its extended error messages. The listing which was 
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produced flagged all areas of the program which required 
revision. Program changes were accomplished utilizing «his 
WATFIV listing. Some of the mcre general and repetitive 


changes are listed in Table III. 


TABLE III 


FLO27 Re-Programing Changes 


CDC Code IBM Code Change 


Mariabiles: ee ee 
FURIT, FWRIF, 
aRear, IWkRIT, IREAD, 


WRITE (IWRIT,600) WRITE (6,600) 
READ(IREAF,500) READ (5,500) 


READ ae WRELE 7 OT Changed to READ 14, 
REWIND 7 WRITE 14 9c REWIND 14 


me(UNIT (N ).GT.0.) GO TC ALL of this type 


Eliminated from program 


iminated 


The most difficult change +o make occurred with the CDC 


Buffer IN or Buffer OUT statements which were used in the 
program to transfer portions of a three-dimensional array 
into and out cf main memory. The Buffer routines reduce th2 
memory size required to execute the program. This statement 
tvyce cccurred inthe main program and severai of the 
subroutines. 
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The change required to translate this statement is presented 
below with the CDC code preceeding the IBM FORTRAN. 

PUFFER, OUT ANS, 1) (6(1,1,1) ¢G(NX,8Y,1) ) changed.to 

WRITE (N3) ((G(I,J,1) ,I=1,MX) ,J=1,MY) and 

EUFFER IN (N1,1) (G(1,1,M),G (MX,MY,M) ) changed +o 

READ(N1,ERR= ) ((G(I,J,M) -L=1,MX) ,J=1,MY) 

The variable ERR was assigned the GO TO statement number 
of the UNIT statement immediately following the BUFFER IN 
fre of code. AS an example, if the UNIT statement 
follewinao the EUFFER IN cede was - If(UNIT(N1).GT.0.) GO TO 
151, then the number 151 was assigned to variable ERR 
follcwing the ¢qual sign. All CDC UNIT statements were 
eliminated from the FLO27 source code per Table II. 

In addition to the pregram changes required to run FLO27 
on the IBM ccmruter, several lines of code were added to 
modify the output format to a more usable form. A 
Subrcutine, VERTEC, which calls the Versatec plotter was 
also Pak d to enhance prcegram usefulness. This plotting 
routine is user controlled through an input variable and is 
ex .lained in the next section. The modified FLO27 progran 


source code is presented in Appendix E. 
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To facilitate program data entry several input variables 
which had reccmmended values were initialized to these 
values within the Main pregram and the subroutine GEOM. The 
initialized input variables and their values are presented 


in Table IV. 


TABLE IV 


Initialized Input Variables 
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A comclete description of each input variable in Table 


IV can be found on pages 19 through 23 of (Ref. 5}. 


B. PROGRAM DESCRIPTION 

The FLC27 program is a computer code written to analyze 
the transcnic flow over awing alcne or a wing ona 
cylindrical fuselage. It uses a finite-volume formulation 
zo sclve the exact potential flow equation in conservative 
form. In the development of the equations, the basic 
assumptions are; eveads flow, no heat or work transfer, 


m@sentropic flow, irrotational flow, no body forces and a 


verfect gas. The velocity vector in cartesian coordinates 
is 
“a a “A 
V= ui + vd + wk (2.1) 
where u, v and w are +te velocity components. The 


continuitv equation, assuming steady flow, is 


o (eu) ee gee (2.2) 
Next a velocity potential is introduced such that the 


velocity components are calculated as <=he gradient of this 


votential. 


u =o pny =o , w =O (2.'3) 
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With the introduction of the velocity potential, the 


continuity equation 2.2 becomes 
XH) * 04) + 34) = (2.4) 


Assuming no heat or work transfer, the energy equation can 


be written as 
Mi/=T {1+ ——— (2.5) 


The flow is assumed to be uniform in the far field. On the 
surface of the bedy, the normal velocity component is zero. 
The velocities and densities of the near field are 
normalized using the free stream velocity and density, thus 
we = 1 and Ceo = Ue Using the assumptions that the flow is 
isentropic and a perfect gas, the energy equation 2.5 can be 


shown to be 


- 1 2 -1 
0 = 41 + oe) M.., (1-4) € (226) 


With equaticns 2.5 and 2.6 there are only two unknowns, 
d and e . They can be solved, subject tc the boundary 
condition of flow tangency, using a finite volume technique. 


The basic numerical scheme for «the solution is the 
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constructicn cf a mesh frem small volume elements (cubes) 
vA 

which are packed around tke wing or wing body configuration. 

The cubes in the computational domain ars separately mapped 


+o distorted cubes in the physical domain by independent 


transformations from local coordinates X, Y* and “Z”° to 
Cartesian coordinates x, y and z. The mesh points are the 
vertices (ccrners) of the mapped cubes. The velocity 


potential and density are calculated at each verctex in the 


mesh. The pressure distribution can then be calculated fron 
ee (2.7) 


In the event that tke local flow velocity becomes 
sucsersonic and shocks occur, these are handled in the usual 
manner by insuring that: 

1) The tangential’velocity components are equal on each 
side of the shock. 

2) Continuity is maintained by keeping the product of 

0 Un constant across the shock (where UJ, is the normal 
velocity component). 

3) Discontinuous expansions (corresponding to an 


“expansicn shock") are excluded from the flow field. 
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The assumoticn cf isentrorfic flow along with the existence 
ef shecks eresents a contradicticn which can only be 
resolved by limiting the flow to very weak shocks for which 
entrepy and verticity generation may be ignored. Thus, 
Solutions will. be valid only for subsonic free stream 
velocities. 

The main three-dimensional array containing the 
potential function data is stored on disk, and special 
unformated incutyoutput statements are used to bring planes 
of data into céntral computer memory and to store updated 
planes of data back On the disks «Inaths construction of the 
computational coordinate system, a Joukowski transformation 
is used to transform the cylindrical fuselage to a vertical 
slit and then a sheared parabolic transformation is used in 
planes containirg the airfoil sections. A detailed 
mathematical formulation of the potential flow analysis is 
contained in [Ref. 7]. 

T. Program Input 

The ingputsto,FELO27 consists of variables wien are 
read with an 8F10.6 FORMAT. Eacheeinput card has a title 
card which precedes it. Thisptitle cand,.contains the ,input 


variakle name and effectively labels the input data for easy 
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reference. The title for each input data is placed in the 


same column as the input data it labels. The*titcie cards 
are read with a 20A4 FORMAT. All numerical input values are 
real numbers. The following data deck, listed card by card, 
is the minimum input data required for a simple wing 
analysis. Each "card" can be interpreted as one line of 
data on your terminal. A complete sample data set is 
presented in Appendix C. 
CARD 1 The Run Title (64 characters maximum) 
CARD 2 Title card for the input variables 
FPNX, FNY, FNZ, FMESH and FPLOT 
CARD 3 
Cols. 1-10 FNX - Number cf computational cells in 
the chordwise direction for the initial mesh. 
MAX = 160/2**n, wkere n = FMESH -1. (See Cols. 31-40 
for FMES&) 
Cols. 11-20 FNY - Number of computational cells in 
the nermal direction £eem the airfoil surface for 
the initial mash. 
MAX = 16/2**n, where n = FMESH -1. 
Cols. 21-30 FNZ - Number of computational cells in 
the scanwise direction for the initial mesh. 


MAX = 32/2*n, where n = FMESH -1. 
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Cols. 21-40 FMESH - Determines the number of times a 


program generated computational mesh is refined. 
ENpLeOs.OnLys ds Up.) 2e0) OF 3-0 fOEYCOaErSe, Medium.or 
fine mesh. If 3.0is selected the pregram will 
calculate flow over the wing for the coarse mesh 
then half the mesh size (medium), recalculate, then 
half the mesh again (fine) and doa final potential 
flow calculation. Output parameters are printed for 
each mesh size fer which calculations were 
performed. 
Cois....41~-50 «FELOT ,=.0utput-flag 
0.0 = Normal cutput without printer-plot of Cp 
1.0 = Normal output with printer-plot of Cp 
at each computational mesh point for 
each wing section. 
2.0 = Normal cutput with Versatec plots of 
Cp versus X/C for each wing section of 
the final mesh. 
Title card for the input variables 
FIT, COVO and P10 
Cne card for each computational mesh. Total 


numker of cards equal to M = FMESH. 
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CARD 6 


CARD 7 


Cols. 1-10 FIT - A parameter which fixes thes 
maximum number of iterations the program will use to 
converge the velccity pctential toa specified 
tolerance (COVO). This parameter must be repeated 
for ¢ach mesh refinement. 

Cols. 11-20 COVO - Velocity potential convergence 
criteria. This input variable is also entered for 
each selected mesh. A value of 0.000001 is 
reccmmended. 

Cols. 21-30 P10 - This parameter determines the 
suksonic point relaxation factor for the specified 
mesh size. A value of less than 2.0 must be entered 
for each designated mesh. Pecommendsd values are: 
126,f¢65 coazse,,. t.3.for wediunm and 1.2. for.the fine 
mesh. 

Title card for the input variables 


FMACH, YA, AL and CDO 


Cols. 1-10 FMACH - Free stream Mach number 


Cols. 11-20 YA - Yaw angle in degrees 


Cols. 21-30 AL - Angle of attack in degrees 
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CARD 8 


CARD 9 


Gols. 31-80 CDO' = Drag coefficient due to skin 


Ericticn. Unless known, an estimated value of 0.01 
is reccmmended. 
Title card for the input variables 


ZSYM, FNS, SWEEP, CIHED and FUS 


Cols. 1-10 ZSYM - The wing planform symmetry 
trigger. 
0.0 = Yawed wing, has no spanwise? symmetry 


1.0 


Swept wing, has spanwise symmetry 

Cols. 11-20 FNS - This input variable <+¢lls the 

oregram the total number of wing sections you have 

selected to define the wing half span. The number 

must ke at least three (3) but not more than eleven 

(11) sections. 

Cols. 21-30 SWEEF - Leading edge sweep angle in 

degrees. 

Cols. 31-40 DIHED - Dihedral angle in degrees. See 
Eide. Zs Vs 

Cols. 41-50 FUS - Input the fuselage radius. Enter 


0.0 for a wing-alone case. 
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Figure 2.1. Dihedral Angle 


Data input cards from 10 through 15 are used for 
defining wirg flanforms and section géeometrics. For the 
feet ae secticn, all data cards from card 10 through card 
15 must be used. For the second and subsequent sections 
there is an cption for Skipping the wing section defining 
Marea (Cards 12 through 15) and copying the data from that of 
the previcus section. This option is controlled by the input 
Variable FSEC. Tf this oftion is not used, data cards from 
10 through 15 must be repeated for each wing defining 
section. The number of wing sections which are defined is 
insut with the variable FNS. Remember, up to 11 sections may 
be defined, anda minimum of 3 sections must be defined. 
All wing planform and section defining geometrics must be in 
consistent units. Wing planform and section defining 


suantities are presented in Fig. 2.2. 
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Figure 2.2. Wing Defining Geometry 


GARD 10 Title card for the input variables 
Sig hie VG, CHORD, THICK, AT and FS EC 
CARD 11 
Cols. 1-10 ZS - The section spanwise coordinate 
(Start at the canterline and work 
outboard) 
Ccls. 11-20 XL - Section leading edge X coordinates 
Cols. 21-30 YL - Section leading edg= Y coordinate 


Cols. 31-40 CHORD - Section chord length 
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CARD 12 


Cols. 41-50 THICK - The thickness scaling factor 


can be used to scale all Y coordinates of the wing 
section. Thus percent thickness and camber are 
increased (or decreased) accordingly. Use: Ont 
no scaling is desired. 

Cols. 51-60 AT - Th? twist angle of each section 
(gecmetric twist) measured from the X axis to the 
chord dine. A positive twist angle reduces the 
section angle of attack and gives "washout". Use 
Q.0 fcr no twist. 

Cols. 61-70 FSEC - This is a flag which determines 
whether or not the program reads wing section 
defining geometry from a previous wing section or 
from new defining geometry. For the first section 
defined you pust set FSEC to 1.0. Following the 
first section, if you define new section geometry 
then use FSEC = 1.0. If you want the program to 
read the section geometry defined from the previous 
secticn then set FSEC = 0.0. 

Title card for the input variable 


FN 
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CARD 13 
Gels. 1-10 FN - This variable contains the number 
ef peints which define the upper and lower surface 
of the section. A maximum of 161 points may be 
used. 

CARD 14 Title card for tke input variables 
LE(L) and YP (2) 

eaepS 15-1 to 15-N Total number of cards equels N, 
where N = integer part of (FN+2)/3. 
The X and Y coordinates at each point are entered 
in pairs, three points to a card. (See Appendix C 
for sampl= input) 
Geiss 1-10 XP (2) = K coordinate of ‘the wing 

section point 
11-20 YP(I) - Y coordinate of the win 


secticn point 


21-30 defining X coordinate for next 
point 

31-40 defining Y coordinate for next 
point 

41-50 defining X coordinate for 


following point 
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51-60 defining Y coordinate for 
following point 

The X and Y coordinates of the wing section defining 
points must be entered starting with the upper surface 
+railing edge peint and preceeding along the upper surface 
to the leading edge, and returning along the lower surface 
to the lower surface trailing edge point. Ee. es Very 
im-ortant to define the section leading edge with a large 
number of closely spaced points. Suggest at least 0.05 
snacing or less ketween X coordinate values from 0.1 X/C to 
the leading edge, X/C = 0.0. Bach Xx and, ¥Y coordinate point 
is nermalized using the chord length for that section. 


Section defining geometrics are illustrated in Fig. 2.3. 
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a chord — 


Figure 2.3. Section Defining Geometr 
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ennp 16 "Title card™-centaining the words in’ Cols. 1-80 


END OF CALCULATICNS 
GARD 17° Title card for the input variable 

FNX 
CARD 18 

Cols. 10 FNX - This variable indicates the end 

of a set of calculations and must be set equal to 

0.0. Its purpose is to indicate that the program 

has run to completion. 

2. Pregram Output 
Outout from the FLC27 program varies with the value 

of the input variable FPLCT. When FPLOT is set equal to 0.0 
a normal output is produced. This normal output contains 
fen crder of occurance) : refined input geometry data 
emcluding trailing edge slope andangle calculations; 
iterative soluticn of the potential flow mesh; section 
characteristics and wing characteristics. The iterative 
Pern Zion, section and wing characteristic data are repeated 
for each mesh refinement requested. Thus, ?+Lti-Dhe+Lnput 
variakle FMESEH is set equal to 3, these data are calculated 
and cutput three times. The last data in the normal output 


consists cf the non-dimensionalized chord (X/C) and pressure 
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coefficient (Cp) data at each computational mesh point for 
each wing section calculated during the final mesh. A 
sample of the normal output data is presented in Appendix D 
and rerresen+s the output data from Appendix C input daté. 

If variable FPLOT is set equal +0 1.0, cthe output 
data is increased considerably. This output contains the 
nozrmal output plus a line printer-plot of theepressure 
coefficient at each computational mesh point for 9ach wing 
section. The line printer-plot is produced for each wing 
section of each mesh refinement. The length of the output 
data with FPLOT set equal to 1.0 can approach 6000 records 
d2cending on the number of mesh refinements requested. 
These vlots are of questionable value and, ‘therefore, an 
alternate plotting program was developed. 

When the variable FPLOT is set equal to 2.0, the 
normal output data is preduced plus a Versatec plotting 
Subrcutine (VERTEC) is called. The subroutine outputs, via 
the Versatec rlotter, plots of Cp versus X/C for each wing 
section of the final mesh calculations. Thas,croutine is 
Simply putting into plot form the Cp and X/C numerical data 
contained in the normal output. A sample of the Versatec 


plot is vresented in Fig. 2.4. 
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Figure 2.4. Versatec Plot of Cp vs. X/C 
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IIT. INTERACTIVE INPUT PROGRAM FLO27IN 


The input data file required for the FLO27 program is 
extensive. Errors in input data FORMAT will cause program 
errors at execution time. In order to eliminate these 
errors and reduce the input data workload, a computer 
F2rminal interactive program was written. Thais interactive 
program, called FLO27IN, is a user-friendly way of creating 
an input data file fcr the potential flow wing analysis 
program FLO27. The FLO27IN program source code is presented 
in Appendix F. 

The interactive program, FLO27IN, when executed ask's 
guesticns of the user in crder tc construct and write to the 
user's "A" disk the required FLO27 input data file. The 
follewing presents the step-by-step procedure for executing 
the interactive program FLO27IN. 

SIEP #1---Log on to any IBM 3033 interactive terminal 
with your user number and password. 
STEP #2---Once logged on and in the CMS operation mode 
type: 
CP LIEK O0247P 191 120 RR then hit ENTER. 
STEP #3---The word PASSWORD will appear, Type and ENTER 


AERO 
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STEP #4---Type and ENTER 


Acc 120 D 
STEP #5---Type and ENTER 
LOAD FLO2Z7IN (START 

The screen will display the header for the interactive 
program. Answer 2ach question presented. At the end of 
each guestion in parenthesis is the input data variable 
associated with that question ard whether the input 
carameter is a real number (R) or an intager (I). Example: 
==> Enter the free stream Mach number (FMACH): (R). FMACH 
is the invut data variable for the question. As you proceed 
through the FLC27IN progran, opportunities to review and 
change input data will be presented. Should it become 
necessary to change your input data after completing the 
FLO2Z7IN srogram, you can simply XEDIT the created data file. 

The FLO27IN program also inccrforates a library which 
contains the wing-section defining data for a number of 
Current wing shaces. A copy of this library is presented in 
Ac-endix B. This feature will be displayed during progran 
execution by the use of a menu from which the user can 


select a pre-defined wing section or define his own. 
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Upon completion of user inputs to the interactive 
vsroaqram three additicnal data lines are automatically 
mratten to'the bettom of the input file. They are: 

END OF CALCULATION 
FNX 
9.0 
anraddition; Job Control Language (JCL) cards are 


meattecen £6 the top and bottem of the file. All JCL cards 


Beadit With a // format. After FLO2Z7IN has run to completion 


ty7e and enter RELEASE 191 to release the aero disk which 


was linked while executing the FLO27IN program. The created 


data fil2 is written to the user's "A" disk with <filename> 


<file+vype> of FLO27 DATAIN. Additional changes can be made 


Simply by entering th2 XEDIT mode and editing the file. 
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Lv. £LO2Z7 BAICH SYSTEM EXECUTION 


The potential flow prcegram FLO27 can be executed after 
the inout data file has been created. The batch processor 
is required fcr FLO27 execution because of the extensive CPU 
time needed to run the progran. While in the XEDIT mode, a 
standard JCB card must be added +o the top of the FLO27 
DATAIN file vorior to submission for job 2xecution. The JOB 
card has the form: 

//jotname JOB (nnnn,pppp) ,tident!,CLASS=J 

jobname = may contain up +o 8 alphanumeric characters, 
the first of which must be alphabetic. 

mnnn = your user number 


pppp = project number, assigned by professor 


fident' = contains the user's own identification 
information. A maximum of 20 characters may be 
contained within th2 single quotation marks. 

After adding the JOB card to your data file, you are 
rsady to execute the program. Type SUBMIT FLO27 DATAIN and 
press ENTER. Eatch runs are normally not worth waiting for. 
To inguire about the status of the job, enter INQ and the 
job name used on the JOB card or "logoff". If the system is 
busy and the maximum mesh size was selected, it may be 


several hours before your job is run. 
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When the jok is run the output will be spooled to «he 
batch printer located next to the VM printer in the main 
eomgruter?P buildings; «The title at the top of the printout for 
batch jobs is the name entered on the JOB card. EE ites 
desired to have the program output data spooled directly to 
the terminal, it will be necessary to add one additional JCL 
card to the input data set. This card must be placed 
immediately following the JOB card and has the form: 

//*MAIN ORG=NPGVM1.nnnnP 
where nnnn = ycur user number 

Inserting this card in the input data will cause all 
program output to be spocled to the user's virtual reader 
where it may ke looked at, printed or transferred to his "aA" 
disk. To enquire as to whether information is in the reader 
Simply type RDR and hit enter, then follow the instructions 


on the screen. 
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V. PROGRAM FEST RESULTS 


The FLO27 program was tested in three stages; (1) during 
the reprograming phase for conversion completeness, (2) 
after sucessful conversicn with suitable wing data for 


program accuracy and (3) during an AE-4501 class project. 


A. ACCEPTANCE TEST DATA 

To test and ensure that “the FLO27 program was converted 
to IEM compatible Fortran without error, an acceptance test 
data set was used. The acceptance test input and output 
data was supplied with the criginal CDC program source code. 
After conversion of +he FLO27 program to Fortran suitable 
for the NES IEM system, the acceptance test input data were 
run and «he output results compared to the output generated 
by the CDC systen. 

Both outout data sets were numerically exact when the 
FLO27 program was run in double precision on the IBM system. 
If the oregram was run ir single precision, the numerical 
output values were exact tc the third décimal place. The 
difference in single precision accuracy occurs because the 
CDC system uses a 64 bit word length while the IBM systen 


mena Jength in single precision is only 32 bits. It was 
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decided thact the IBM single precision accuracy was 


Satisfactcry. 


B. CCMPARISCN WITH OTHER FROGRAMS 

The FLC27 program was also tested for accuracy by using 
the wing vlanfcrm and section data from a NACA 572 wing. 
The data were run on both the FLO27 program and the Douglas 
potential flow program [Ref. 1]. The data generated by both 
programs was compared to wind tunnel data for the NACA 572 
wing (Ref. 8]. The results are presented in Fig. 5.1 as 
prots cf lift ccefficient versus angle-of-attack. The 
results show that for the NACA 572 wing the FLO27 program 
more accurately predicts the wing lift coefficient than does 


the [Couglas program. 


Cc. AE-4501 CLASS PROJECT 

The final test phase was conducted by introducing the 
FLO27 program into the AE-4501 course as a class project. 
This was accomplished to determine student problems/comments 
concerning the data input program FLO27IN and to test an 
additional wing shape. The wing chosen for study was that 
of the A-7 airplane. The A-7 wing has a distinct leading 
edge netch at the approximate mid-span. When the planform 


geometry was run with the notch included the FLO27 program 
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Figure 5.1. Program Calculated and Wind Tunnel Data 


Tan to completion but gave negative values for section and 


total induced drag coefficient. The value for the lift 
coefficient was low for the freestream Mach and 


angle-of-attack used. Et Was found that if the notch was 
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excluded from the wing gecmetry input data the program 


results were satisfactory both for induced drag and lift 
coefficient. 

From the AE-4501 class experience it was determined «hat 
sharp wing clanform discontinuities cannot be handled by the 
program. If hcwever, the changes in shape are gradual, such 
as awing glove, the program output appears to be 
Satisfactory. Such was the case with the acceptance test 
case data where the wing geometry was that of the F-8 


supercritical wing which incorporates a wing glove. 
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APPENDIX B 


LIERARY OF AIRFCIL SECTION GEOMETRIES 


user input section coordinate data 


flat 


plate data 


symmetrical wing (11% thickness at 30% chord) 
12% thickness at 32% 


supercritical wing (cambered, 
chord) 


NACA 24-30-0 (cambered, 
F-14 wing (cambered, 


12% thickness at 30% chord) 
9.5% thickness at 37% chord) 


A-7 wing (7 deg droop at 20% chord, 7% thickness at 
43% chord) 


LISSAMAN 7769 Airfoil 


chord) 


NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 


0010 


0010-35 


0010-64 
0010-66 


16-009 
63-010 
63A010 
64-010 
644010 
65-010 
65A0 10 
66-010 


(Symmetrical, 
0010-34 


(symmetrical, 
(symmetrical, 
(symmetrical, 
(Symmetrical, 
(Symmetrical, 
(Symmetrical, 
(symmetrical, 
(Symmetrical, 
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